Quantum theory of reactive collisions for potentials 
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We develop a general quantum theory for reactive collisions involving power-law potentials 
(— valid from the ultracold up to the high-temperature limit. Our quantum defect framework 
extends the conventional capture models to include the non-universal case when the short-range 
reaction probability P^^ < 1. We present explicit analytical formulas as well as numerical studies 
for the van der Waals (n = 6) and polarization (n = 4) potentials, and relate our work to recent 
merged beam experiments on Penning ionization for collision energies below 0.1 meV [Henson et al. 
Science 338, 234(2012)]. 
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Much recent work involves the inelastic and reactive 
collisions of cold atoms or molecules with one another [Ti- 
ll] or with ions in hybrid traps [5]-[7] . These could involve 
the ultracold regime with translational temperature on 
the order of /iK or less or the cold regime between a 
few mK and a few K. Systematic theoretical principals 
for understanding the quantum dynamics of such colli- 
sions are needed. Much work has already been done in 
this area 5 as reviewed by Ref. [8 . One class of theories 
based on Quantum Defect Theory (QDT) allow us to sys- 
tematize and develop tools for understanding such colli- 
sions [9 -14 . One special limiting case is that of highly 
reactive collisions, where simple classical trajectory cap- 
ture models known as the Langevin (n = 4) [15] or Gorin 
(n = 6) [16 models apply when the long range potential 
takes on the form —Cn/R^ (n > 3). These familiar mod- 
els assume that every classical trajectory contributes to 
the collision cross section that is captured by the long 
range potential so the particles spiral in to short dis- 
tance where they react or relax with probability P^® . We 
use the term " universal" to describe capture models with 
P^® = 1, since they do not depend on any details of the 
strong short-range chemical interactions. 

In the cold and ultracold regimes, it is essential to build 
in quantum corrections to these classical models due to 
quantum threshold laws [ TTHTO] . This has been done us- 
ing QDT for both the Langevin [20l|2T] and Gorin [221 [23] 
universal models where P^® = 1. We will follow here the 
QDT formalism by Idziaszek et al. [211 HI] which also 
treats the nonuniversal case where P^® < 1 and gives an- 
alytic formulas for the threshold properties as T 0. 
This theory in its numerical form allows us to calculate 
and interpret the ionization rate constants measured by 
recent merged beam experiments in the cold regime [24 . 

Generalized complex scattering length. We consider 
reactive collisions of particles interacting via power-law 
potential V{r) = —Cu/t^ {n > 3) at large distances 
r > Po, where Rq denotes the range of short-range 
forces. These include, for example, collisions of 5'-state 
atoms or molecules in the rovibrational and electronic 



ground state (n = 6), or collisions of 5'-state atoms 
with ions (n = 4) [25H27j . The characteristic length 
Rn = (2/ia/n2)'/^""'^ and energy = h^/(2^Rl) 
are associated with the long-range potential. Typically, 
En is of the order of the centrifugal barrier height for p- 
waves, while R^ defines the characteristic distance for the 
position of the barrier. The scattering channels are de- 
fined in terms of their internal states and the partial- wave 
quantum numbers Im. For simplicity we will use only 
Im quantum numbers to specify the channels, while the 
internal quantum numbers will affect only the quantum- 
statistical prefactor which will be discussed later. Fol- 
lowing [22j ^ we define an energy-dependent [29l [30] 
complex scattering length for each m 
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defined in terms of the diagonal elements of S matrix. 
The elastic /C®^ and the reactive JC^^ rate constants can 
be expressed as ^ 

IC^\E) = Y.Kf^{E) = 5^ ^ |1 - S,^,,m{E)\' , (2) 



1C^%E) = ^KUE) = 5- ^ (1 - \SemME)n , 



(3) 



where k'^ = 2iiE /h? with E denoting the total energy, /i 
the reduced mass, and g is equal to 1, except the case 
when the particles are indistinguishable and are in iden- 
tical internal states. In the latter case, g = 2 and ^ is 
restricted to even (odd) numbers for bosons (fermions). 
Comparison with experimental data requires averaging 
over actual distribution of energies in the particular ex- 
perimental setup. 

Quantum Defect Theory (QDT). We adopt the formu- 
lation of QDT by Mies [TOl |TT] to find a general formula 
for aim- Mies uses different QDT functions than in Gao's 
similar treatment [14 , thus allowing additional insight 
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into the dynamics. The basic idea of QDT is to sepa- 
rate short-range and long-range properties of the wave 
function. To this end one introduces a pair of hnearly 
independent solutions fir^E) and gir^E) for each par- 
tial wave I that have local WKB-like normalization at 
short distances, as done in [22]. One then parametrizes 
the wave function at short range using / and g functions 



^{r,E)=A{E) f{r,E) + Y{E,£)g{r,E) 



(4) 



Here, A{E) is the amplitude, and Y{E, i) is the quantum- 
defect parameter, that plays a key role in the quantum- 
defect method, parametrizing the short-range behavior of 
^(r). In many cases when short-range forces act at dis- 
tances much smaller than i?^, one can neglect the energy 
and angular- momentum dependence of Y{E^£) [2T| [3T]. 
Then we can represent Y in the form Y = —iy^ where y 
is a parameter related to the probability of reaction [22 . 
Using a WKB-like representation of / and g yields 
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where k{r) = {E — Ui{r))/h is the local wave vec- 
tor, Ui{r) = V{r) + h^£{£ + l)/{2/ir^), and (3{r) = 
J^dxk{x) is the WKB phase. The first term repre- 
sents the incident flux of particles, and the second one 
the reflected flux. We notice that y = 1 corresponds 
to unit reaction probability (no outgoing flux) while for 
y = there are no losses (incident and reflected fluxes 
have equal amplitudes) [22]. The reaction probability is 
= 1 - [(1 - y)/(l + y)f = 4y/il + yf. 
Using MQDT formulas connecting S with Y{E^ i) p!Ol, 
we derive an exact expression for the energy-dependent 
complex scattering length, which is valid for arbitrary 
energy and partial wave i 



a£m{E) = - ^ tan 
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Here, £^{E^i) is the scattering phase shift of the poten- 
tial V(r), and C{E) and tanA(£;) are the MQDT func- 
tions that connect the solutions with short- and long- 
range normalization: /(r) = C~^{E)f{r) and g{r) = 
C{E)[g{r)^tdinX{E)f{r)] (see Refs. ^ [ni[22] for de- 
tails). Substituting ([6| into (|3| we can express the reac- 
tion rate directly in terms of MQDT functions 
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(1 + yC-^ {E, £)y + y^ tan^ X{E, i) 
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The C ^ and tan A functions build in the effect of shape 
resonances due to quasi-bound states inside the long- 
range centrifugal barrier Ui{r). References [11] [211 [32 
give examples of these threshold functions. 



It is interesting to investigate the low- and the high- 
energy limits of the total reaction rate = ^g^^Ym- 
At low energies in case of distinguishable-particle col- 
lisions or indistinguishable boson-boson collisions, the 
main contribution comes from 5-wave scattering, while 
for identical fermions scattering is dominated by the p- 
wave. We introduce 
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which is the mean scattering length for — 1/r" potential 
[33] , and the characteristic p-wave volume 



V 



TT (n - 2)("-8)/("-2) 



r2 



(A) 



(9) 



The key result of this letter are the respective s and p 
reactive rate constants in the threshold regime, given by 
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where u = cot and s = a/a with a denoting the 5- 
wave scattering length. One can show that for a fixed re- 
action parameter y, the maximal reaction rate is obtained 
for s ±oo {£ — 0) and for s cos {£ = 1), which 
correspond to the resonant case of a bound state crossing 
the threshold. In the universal regime (y=l) Eqs. (10)- 
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(11) reduce to JCqq 
independent of s. 

At high energies where many partial waves contribute 
we first derive an approximate expression correspond- 
ing to the classical limit of the scattering. We assume 
C~'^{E^i) = 1 and tanA(£^,^) = for partial waves at 
which the collision takes place above the barrier, while 
for collisions below the barrier we take C~'^{E^£) = 0. 
This neglects the effects of the quantum tunneling and of 
the quantum reflection. In this approximation we obtain 
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where ^max(^) is the maximal angular momentum at 
which the top of the barrier is equal to the collision en- 
ergy E. For a power-law potential V{r) = —Cn/r^ this 
leads to 
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This can be reproduced by solving the classical problem 
of scattering on — Cn/r^ potential, assuming that all tra- 
jectories that fall into the collision center contribute to 



the total reaction cross-section cr^® with the probabihty of 
reaction P^®. Then, /C^® = cr^®'^, with v = hk/ ji denoting 
the mean relative velocity in the gas. 

A more accurate treatment requires inclusion of the 
quantum effects, in particular when the collision en- 
ergy is close to the top of the barrier. In our ap- 
proach we approximate the top of the centrifugal bar- 
rier with a parabolic barrier with imaginary frequency 
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jl^En + 1)] '"-^ . For such a potential 
one can find an analytic solution in terms of parabolic 
cylinder functions [34 . At large distances from the 
parabolic barrier one can retrieve the WKB-like form 
of the wave function and identify the components cor- 
responding to incoming and outgoing fiuxes of particles. 
Parameterizing the wave function at small distances in 
terms of the y and s parameters (cf. Eq. ([5|), we calcu- 
late the S matrix and the loss probability 1 — \S£rn{E)\'^. 
As we are interested only in the behavior of thermally- 
averaged reaction rates at large energies where several 
partial waves contribute, we average the result over the 
short-range phase. For a parabolic potential we derive 
the following formula, valid for arbitrary energy 



P,,{i,y) = l-\S,m{E)\^ = 



(14) 



where e{i,E) = E/huji - V2n - 4.^^^^^^ is the dimen- 
sionless energy measured with respect to the peak of the 
parabola in harmonic oscillator units. In the universal 
regime y = I Eq. ( 14 ) is identical to a textbook solution 
derived in | 

by 



The reaction rate at high energies is given 
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where Yl' denotes summation over angular momenta al- 
lowed by symmetry. At high energies we can replace 
summation over i by integration, which yields Eq. (16). 



The quantum-statistical factor g disappears since for in- 
distinguishable particles the summation is done only over 
even or odd values of £. 

The loss probability P\s{i,y) as a function of a con- 
tinuous variable i{i -\- 1) is shown in Fig. [l] The figure 
compares the loss probability calculated in the parabolic 
potential approximation with the classical approach as- 
suming that only collisions with energies above the bar- 
rier contribute to the reaction rate. The latter exhibits 
a step-like behavior, while the former resembles a Fermi- 
Dirac function. We observe that for y = 1 classical de- 
scription overestimates the reaction rate in the regime 
affected by the quantum refiection (energies above the 
barrier, marked in orange), and at the same time does 
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FIG. 1: Reaction probability calculated for a parabolic po- 
tential fitted to the actual centrifugal barrier of the van der 
Waals potential as a function of a continuous variable 1) 
(red solid and green dashed lines). The case P^^ is averaged 
over the short-range phase and is thus independent of s. The 
Langevin approximation depicted by black solid lines assumes 
constant reaction probability P^^ above the barrier, and no 
reaction below the barrier. 



not include the contribution from the quantum tunnel- 
ing (energies below the barrier, marked in blue). In the 
universal regime y = 1 both contributions turn out to 
be almost equal, and in this particular case the classi- 
cal description should work relatively well. In contrast, 
for y < 1 the loss probability is additionally affected by 
the shape resonances which makes the contribution from 
the quantum tunneling typically larger. In such a case 
the two effects do not cancel each other and the classical 
theory underestimates the reaction rate. Typically, for 
weakly reactive particles the quantum correction with re- 
spect to the classical theory is of the order of few percent 
at E^ IO^Eq for y < 0.01. 

Examples for van der Waals and the polarization po- 
tentials. Fig. [2] shows reactive rates for van der Waals 
potential for bosons and fermions for two exemplary 
values of s in the universal {y = 1, considered also 
in [20I |23]) and non- universal {y = 0.1 and y = 0.01) 
regimes. The numerical results (bold solid and dashed 
lines) were obtained by imposing QDT boundary con- 
ditions specified by Eq. ([5| at small distances and then 
propagating the wave function to large distances. Af- 
ter extracting the S matrix we calculate the reactive 
rates from Eq. (|3| and then perform thermal averaging: 
{Kre)thiT) = 2/V^ikBT)(^/^^ J dEVEe-^^'"^^Kre{E). 
We note that approximate result (15) works relatively 
well in all cases for energies E > lOO^^e, except for 
y = 1 when there is no contribution from the shape reso- 



nances. Moreover, we have verified that at large energies 
the rates for bosons, fermions, and distinguishable par- 
ticles become equal which agrees with the approximate 
result of Eq. (15). Figure [s] shows the reactive rates for 
distinguishable-particles collisions in the polarization po- 
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FIG. 4: Measured reactive rate constants (black points), 
QDT fit (blue straight line), contribution from ^ = 5 (dotted 
line) and other partial waves (dashed line) for the collision 
of Ar with He*. Experimental data is taken from [24]. The 
fit parameters are y — 0.007 and s = 3. For this system 
Rq ^ 40 ao (ao is the Bohr radius) and Ee/ks ~ lAmK. 



FIG. 2: Thermally averaged reaction rates for van der Waals 
potential, calculated for different reaction amplitudes y and 
short-range s parameters, for fermionic and bosonic particles. 
The dotted lines depict the high-energy approximation (13) 
assumed in the Langevin theory. Th e dot-dashed lines show 
the high-energy approximation (15) including the quantum 



effects derived by approximating the centrifugal barrier by a 
parabolic potential. The universal y = 1 curve agrees with 
the result of l23l. 
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FIG. 3: Thermally averaged reactive rate constants for atom- 
ion potential, calculated for different reaction amplitudes y 
and short-range s parameters, along with the Langevin ap- 
proximation (dotted line) and our high energy treatment (dot- 
dashed lines). 



tential, where the high temperature limit is described by 
the Langevin theory [151 ES] • For n ■ 



4 Eq. ([13|) predicts 
a constant reactive rate. At energies E > lOO^^e we find 
a good agreement of our numerical results with the ap- 



proximate result (15) based on the parabolic potential. 



Penning ionization. Finally, we apply our theory to 
the case of Penning ionization of argon by metastable 
helium, measured using merged molecular beams in the 
cold regime [24]. The long range part of the Ar-He"^ in- 
teraction is given by van der Waals potential. The short 
range part can give rise to energy- and ^-dependent cor- 
rections to y and s. However, we found good agreement 
with the experimental data at energies below ~ 100 
using the simplest form of our model using the same y 
and s for all partial waves. The position of the reso- 
nance is mostly determined by the s parameter, while y 
sets its width and height. As a result of fitting we ob- 
tained y = 0.007 and 5 = 3, and the scattering length 
a = (1.43 — 0.017i)RQ. The experimental data and the 
theoretical model are compared in Fig. [4j The calcula- 
tions shown in Fig. |4] were averaged over the velocity dis- 
tribution of a merged beam. We found that the observed 
resonance is due to the contribution from the single par- 
tial wave £ = 5. 

In conclusion, we present a simple yet general model 
based on QDT of a reactive collision between two parti- 
cles interacting by a power-law potential. We have shown 
how the low energy collisions are affected by shape reso- 
nances and how they contribute to the high energy col- 
lisions. Our theory is the first approach to nonuniversal 
collisions at low energies which does not require calcula- 
tions of potential energy surfaces and can be applied to 
various kinds of problems, such as molecular collisions, 
atom-ion collisions and Penning ionization. 
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